Entropic uncertainty relations for the ground state of a coupled system 
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There is a renewed interest in the uncertainty principle, reformulated from the information the- 
oretic point of view, called the entropic uncertainty relations. They have been studied for various 
integrable systems as a function of their quantum numbers. In this work, focussing on the ground 
state of a nonlinear, coupled Hamiltonian system, we show that approximate eigenstates can be 
can be constructed within the framework of adiabatic theory. Using the adiabatic eigenstates, we 
estimate the information entropies and their sum as a function of the the nonlinearity parameter. 
We also briefly look at the information entropies for the highly excited states in the system. 

PACS numbers: 05.45.Ac, 03.67.-a, 05.30.-d 



I. INTRODUCTION 



The uncertainty relations, that express the inability to 
simultaneously measure the states of two non-commuting 
observables, form the cornerstone of quantum physics. 
For any pair of operators A and B, in standard form 0, 
they are stated as, 



AA AB > -(■9\[A,B]\^) 



(1) 



where and AB represent the dispersions in A and 
B and [A, B] is the commutator. In recent years, there 
has been a revival of interest in the uncertainty relations 
reformulated from the stand point of information theory, 
called the entropic uncertainty relations (EUR) 0]- For 
instance, the position-momentum uncertainty relation is 
formulated as follows; Given an eigenstate of a quantum 
system, ipil) s-nd tp{p), in configuration and momentum 
space representations, and if Sq and Sp represent their 
information entropies, then the entropic uncertainty re- 
lations can be written down as. 



Sa 



Sp ^ ^QP 



(2) 



where Sqp is the lower bound to the entropic sum or EUR. 
Here, the information entropy is defined as, 



(3) 



where p{^) = is the probability density. The in- 

formation entropy is a measure of the spreading or local- 
isation of the given eigenstate. Apart from its intrinsic 
value, the reformulation also seeks to address some of 
the shortcomings in the standard statement of the un- 
certainty principle ^2]. It also quantifies the uncertainty 
more accurately than the standard statement based on 
dispersions 3]. In particular, a good amount of work has 
focussed on obtaining the lower bounds Sqp in general 
and we mention the result due to Bialynicki-Birula and 
Mycielski 3, 



where D is the dimensionality of the system. Inspired 
by this general result, several authors have focussed on 
obtaining the lower bounds for EUR of quantum systems 
whose classical limit is integrable, like the particle in an 
infinite well |5|], harmonic oscillator and hydrogen atom 
in one and higher dimensions Q , power-law wavepackets 
and oscillating circular membrane @. Apart from 
the lower bounds, the values of EUR as a function of 
quantum numbers has been analysed in detail in these 
series of papers. Recently, Dehesa et. al. have shown 
that for the one-dimensional power-law potentials of the 
form V{x) = x'^^, where is a positive integer, in the 
region of highly excited states, the entropic sum goes as 
ln(2n) for all A; 0. 

On the other hand, the single-particle probability den- 
sity p{^) is also the quantity of fundamental interest in 
the density functional theory and hence its characterisa- 
tion using information entropy as a measure for spreading 
has assumed special interest. In fact, treating atomic 
and molecular entropic sum was considered by Gadre 
[To| and he derived an approximate expression for the 
entropic sum within the Thomas-Fermi framework for 
neutral atoms. Now it is known from several empirical 
studies for atomic, molecular and nuclear distributions 
that the entropic sum can be modelled as JJj, 



S=ci+C2 IniV 



(5) 



Sq + Sp> 15(1 + In vr) 



(4) 



where ci and C2 are constants and TV is number of elec- 
trons or nucleons, as the case maybe. The functional 
form given above seems to be fairly universal for many- 
fermions in some mean interactions 

Thus, one branch of the work on EUR has focussed on 
quantum systems in the classically integrable limit whose 
eigenstates are analytically known. The other comple- 
mentary branch has explored the complex atomic and 
molecular systems using a combination of approximate 
analytical and empirical methods. In these cases, the 
EURs have been obtained as a function of increasing 
quantum numbers. However, simple and chaotic model 
systems that bridge this divide between the purely inte- 
grable and the complex many-body systems have not yet 
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FIG. 1: The contours of the potential as a function of x and 
y for a = 0.2 and k-i,k2 — 1. Each contour represents a 
particular value of energy. For large values of energy, the 
potential develops 'channels', as seen in the figure. 



chaos since they qualitatively mimic the complex dy- 
namics in problems like the hydrogen atom in strong 
fields and generalised van der Waals potential 
a generic chaotic system, most of its eigenstates could 
be thought of as random waves except for a small sub- 
set of 'regular' states that are influenced by the nature 
of local classical dynamics. Hence in this work, we will 
only focus attention, for most part, on the ground state 
and, briefly, on the regular states. We stress that even 
though the system is chaotic, we study only those eigen- 
states that are associated locally with regular classical 
structures. Ground state is pre-eminent because previ- 
ous studies have shown that the ground state saturates 
the EUR inequality in some systems (e.g. harmonic os- 
cillator) and often in complex systems, the focus is on 
the ground state properties. Next, we obtain the ground 
state entropies, as a function of a. 



been considered. The main reason seems to be that, as 
yet, no straightforward analytical technique is available 
to determine their eigenstates. Thus, the main purpose 
of this article is to show that using adiabatic technique 
approximate eigenstates can be constructed and further 
that the entropic sums can also be estimated. In contrast 
to earlier works, we study the EUR as a function of the 
nonlinearity parameter in the system. Thus, we hope to 
understand the effect of nonlinearity using information 
theoretic measures. 

In general, the non-integrable Hamiltonian systems 
are generic rather than an exception and model realis- 
tic physical systems. Their phase space presents a mix- 
ture of regular and irregular trajectories, whose proper- 
ties change with a parameter. The evolution of the prob- 
ability density as a function of the parameter reflects the 
influence of these classical objects. Thus it is important 
to look at the EUR upon variation of the nonlinearity 
parameter. Note that scale invariance of entropic sums 
leads to EURs that are independent of scaling param- 
eters in an integrable system. But, in a non-integrable 
system, the entropic sums could depend on the nonlin- 
earity parameter. 



II. MODEL HAMILTONIAN 

In this work, we will consider the model Hamiltonian, 

^^=f + f + ^ + ^ + «-V (6) 

whose potential is displayed in Fig. ^s-ud ki,k2,a are the 
parameters. The system is integrable for a = and cor- 
responds to two-dimensional harmonic oscillator whose 
asymptotic (large quantum number) entropic sum was re- 
cently derived jQ]. For a choice of a = 0.05 (fci — k2 ~ 1), 
the system displays classical chaos for energy i? > 15 
[T^ . The coupled oscillators are popular as models of 



III. ADIABATIC THEORY 

In the context of two degrees of freedom systems, if 
the frequency of oscillations between both the degrees 
of freedom differ vastly, then such a classical system be- 
comes an ideal candidate for the adiabatic treatment. 
This is the well known Born-Oppcnheimer approxima- 
tion in atomic physics 14\. The approach is to average 
over the 'faster' degree of freedom and incorporate it into 
the Hamiltonian for the slower degree of freedom. In this 
adiabatic approach, the Hamiltonian is effectively inte- 
grable and can be quantised. This has been successfully 
implemented earlier to estimate the energies of the lo- 
calised states in coupled oscillator systems ^2 ■ -^'o'" 
the Hamiltonian in Eq. Certain and Moiseyev con- 
struct adiabatic eigenstates formally Prange et. al. 
construct adiabatic eigenstates for a class of chaotic bil- 
liards [13. 

A. Position Eigenstates 

In this section, we will obtain the ground state of the 
Hamiltonian in Eq. (|SJ), for a << 1, in the adiabatic 
approximation, closely following the technique in Refs. 

We assume that the y- motion is faster and aver- 
age over the fast motion by dropping terms that depend 
purely on x and Px- This gives, 

H'ix;y)^^+'^+axV^E'{y) (7) 

where the variable x is treated as a parameter. This is 
just the harmonic oscillator whose frequency is 

Qx = ^kl + 2ax^ (8) 

and for a « 1, such that fc| >> 2ax^ , then 

VLx'^ uJx ^ k2 + j-x^ (9) 
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The ground state of H'{x\ y) for a << 1 is, 

^{x; y) = (^) ^'^ exp(-c^,yV2) (10) 
Now, the adiabatic Hamiltonian is obtained as, 



Had — — H r LOxJi 



2 2 



x'Jy 



(11) 



where Jy is the classical action corresponding to the 
faster motion. In terms of action-angle coordinates, we 
get, 



Had — Jx 



kl + 2aJy/k2 + k2Jy (12) 



Note that semiclassical quantisation of Had, which is also 
exact in this case, gives for the ground state energy, 



Ead = ^Jkf 



2 



(13) 



The energy estimated by this formula is in good agree- 
ment with the computed energies for small a. Since we 
are interested only in the eigenstate, the last term in Eq. 
(|12|l can be omitted since it only serves to shift the energy 
scale. The ground state of Had is. 



tp{x) = (^-j exp(-cja;V2) 



(14) 



where uj — ^kf + The ground state of the Hamilto- 
nian in Eq. lO, for a << 1, in the adiabatic approxima- 
tion, can be written down in its characteristic form as, 
"^{x,y) = 4>{x]y) ■0(2;), so that. 



(15) 



At this point, in order to be able to do the integration, 
we approximate the integrand for a « 1. We take. 



exp 



-pI\ /-p'(i-(«/fci))a^' 



2ujx 



exp 



2fc2 



k2 



(18) 



Then, the integral reduces to a Fourier transform of a 
Gaussian, 



V2n " \ 2k 



> ' fix) e^"-" dx (19) 



where, 

f{x) = exp 



-x' 



ujp = (3 ■ 



^2 



(20) 



and P ^ ui + a/2k^. This integral can be performed and 
we obtain the momentum eigenstate as. 



1 



^/^^ \k2 



V4exp(-i^)exp(-|L) 



/Wo 



(21) 



This eigenstate is also of the product form charac- 
teristic of the adiabatic approximation, ^(pxTPy) = 
Xi{px]Py)X2{Py), where xi and X2 represent the eigen- 
states of fast and slow degrees of freedom. In contrast 
to the position eigenstate in Eq. H15(l . the ground state 
in momentum representation is normalised only to 0{a). 
This is due to the approximations that were done in the 
course of obtaining the Fourier transform. As a limiting 
case, for a = 0, it reproduces the 2D momentum eigen- 
state. 



Note that '^{x,y) is a coupled eigenstate; the couphng 
provided by lUx- It is correctly normalised. For a = it 
reduces to a two-dimensional harmonic oscillator wave- 
function. The price we pay in the adiabatic approxima- 
tion is that the eigenstate is not symmetric under a; <-> y 
as required by the potential. 



B. Momentum eigenstate 

We calculate the eigenstate in the momentum repre- 
sentation by taking Fourier transform of the eigenstate 
given in Eq. (|15|l . We perform the following integral. 



^'(Px,Pa) = ^ 

Zn 



OO /'CO 



■i>{x,y)e'P-'=e'Pyy dx dy (16) 



— OO — OO 



By a straightforward transformation of variables, we can 
reduce it to an one-dimensional integral, 

"^{Px^Py) = / dx (17) 



%/27r 



1/4 



IV. GROUND STATE ENTROPY 

A. Position Entropy 

Now it is possible to evaluate the information entropies 
with the eigenstates in Eqs. H15I21(I . The position den- 
sity, from Ea. H15|l . is given by. 



p(x,y) = \^{x,y)\^ 



(22) 



Then, the information entropy can be evaluated as, 

poo pOO 

Sg{a,ki,k2) = - / / p{x,y)\np{x,y) dx dy 



— OO J —OO 



h 



(23) 



1 + In TT — In \/uj - 



where Ii is the integral given by. 



/I = / e ^ In ( fc2 -I- ] dx (24) 
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Though this integral can be exactly done, we expand the 
integrand to first order in a and arrive at the following 
expression for the entropy valid for a << 1, 



2.3 



Sq{a, fcl, fc2) = 1 + InTT 



(25) 

The spatial entropy decreases linearly with a. Setting 
a = 0, we recover the correct 2-diniensional harmonic 
oscillator entropy. In Fig Elb) we show the spatial en- 
tropy (fcl = fc2 = 1) as a function of a. For comparison, 
we numerically estimate the entropies by diagonalising 
the Hamiltonian in Eq. © in harmonic oscillator basis 
set and a subsequent entropy calculation. There is a good 
agreement between the numerical results (dots) obtained 
without any approximation and the analytical formula 
(solid line). A linear fit to the numerical entropies gives 
a slope of -0.4533, not far from the theoretical result -0.5. 
The fitted slope closely approximates -0.5 as a ^ 0. 



B. Momentum Entropy 

The momentum density, 'p{px,Py), from Ea. H21|l . is 
given by, 



_ 1 r—r exp 



(-g)»p(-S) 



OJv 



(26) 



Once again, an entropy integral similar to Ea. H23|) can 
be performed by Taylor expanding the denominator in 
the integrand to 0{a). The final, rather cumbersome, 
expression turns out to be, 



Sp{a, fcl, fc2) = 



(1-c) 



1 



-In 



TTf3 



(1 - 2c) 



2(1 -c) 
(27)' 

where we use the shorthand, c = a/2/3fc|. As a limiting 
case, we get for a = 0, the correct momentum entropy of 
2D harmonic oscillator. To obtain an explicit expression, 
we expand to 0{a) all the terms in the above equation. 
The result valid for a << 1, is. 



Sp{a, ki,k2) — I + In tt -|- In 



a fcl 
4 



fc2 



(28) 

This expression shows that the momentum entropy in- 
creases linearly with a. This is evident from Figl^a) in 
which a good agreement is seen between the numerical 
and the analytical result. A linear fit gives the slope to 
be 0.4765, close to the expected value 0.5. 



C. Entropic sum 



Prom the entropy expressions in Eqs. (|25I28|I . we ob- 
tain the sum of position and momentum entropies as, 



(a) 



(b) 
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FIG. 2: Entropy for the ground state of the Hamiltonian in 
Eq. JSJ with fcl = A;2 = 1 as a function of a. (a) Momentum 
entropy and (b) spatial entropy. Dots are numerical values of 
information entropies obtained without any approximation. 
Continuous lines are theoretical curves of Eqs. H25I28^ 
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FIG. 3: Entropy sum for the ground state of the Hamiltonian 
in Eq. Q with fci = fc2 = 1 as a function of a. Dots are 
numerical values and continuous line is the quadratic fit to 
the numerical data. 



S{a)^Sq + Sp^2{l + \mT)+0(a^) (29) 

We point some salient features of this result. To 0{a), 
the entropic sum is invariant under change of a, but def- 
initely satisfies the uncertainty inequality in Eq. Q). As 
expected, due to the scale invariance, it is also indepen- 
dent of the parameters fci and fc2. However, as revealed 
by numerical computations in FigO the significant con- 
tribution to entropic sum comes from terms of 0{oP'), 
which could not be determined from the present method. 
The important result is that, for a << 1, the entropies 
and their sum depend only on the parameter a, infact 
on terms of 0(0?) or higher. Since our method does not 
yield information on the quadratic and higher order a 
terms, we perform a fit to the numerical entropic sum. 
The equation of best fit, for the range of a considered, 
gives S = 4.289 -I- 0.063a -f QAQ&a^. As expected, the 
fitted equation shows that the linear term in a carries 
less weight and the intercept is a good approximation to 
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2(1 + ln7r). The ground state saturates the entropic sum 
only for a = 0. 

The role of a, primarily, is to change the shape of the 
potential which leads to other qualitative changes in the 
classical dynamics of the system. As a increases, the po- 
tential develops 'channels' (see Fig. ^ but the ground 
state does not occupy the channels. The potential re- 
stricts the ground state to remain within the central re- 
gion and not expand into the channels of the potential. 
Hence the spatial entropy decreases with increase in a 
indicating stronger localisation of the 'particle' near the 
origin of the potential. Correspondingly, the momentum 
entropy increases. But they do not cancel one another 
exactly, and the increase in entropic sum is of the order 
0{a^). Thus, it is a rather slow but definite increase, as 
evident from Fig. O 



LARGE PARAMETER 



where U{, , ) is the confluent Hypergeometric function. 
Now, we optimise Ey to obtain b. Hence, 



dEy 



1 kl e'''=2/4" kl Koibk^/Ac 
4 ^ 462 ^ 



4v267ra 
[/(-l/2,0,6fc^/2a) 



(35) 



where Ko{.) is the modified Bessel function of second 
kind. 

Now, we solve for b by setting dEy/db = 0. Once again, 
solving for 6 in a general case does not seem feasible. 
However, in the limit a » 1, we use the asymptotic 
forms as z = (6fc|/2a) ^ [3, 



1 



27- 1 



Koiz/2) 



-ln(z/2) 



(36) 



We note that the limit a — > c» is a singular limit for 
the Hamiltonian in Eq. 0. It does not seem feasible to 
extend the present method to this limit. Hence we look 
at a » 1, but not a — *■ 00. 



A. Position Eigenstate 



where 7 = 0.5772... is the Euler's constant and ^'(.) is the 
digamma function. We use these asymptotic forms and 
ignore terms that contain a in the denominator. Then, 
we obtain. 



dEy 



1 

4 



1± 

Ab^ 



2V2^63/2 



(37) 



For a » 1, we cannot Taylor expand as done in Eq. 
0. Hence, now the adiabatic Hamiltonian will be. 



H, 



ad 



2^2 



kfx 



(30) 



If the terms are rearranged, this turns out to be a quartic 
algebraic equation in b. 



2klb' ^^ + kt 

TT 







(38) 



Presently, exact solutions are not known for this problem. 
To go further, we quantise the y degree of freedom exactly 
and treat the x degree of freedom using the variational 
method. This leads to, 



pI kjx^ 1 



^\lkl + 2aa;2 



(31) 



where Vl^ is substituted from Eq.JHJ. Now, we will use 
variational method to obtain the ground state. The trial 
wavefunction of the form. 



1/4 



exp(-6x2/2) 



(32) 



is assumed. Note that the trial wavefunction is similar 
to the adiabatic ground state in Eq. H14|) . except for oj 
being replaced by the variational parameter 6, to be de- 
termined. Since we are applying the standard variational 
technique as expounded in many quantum physics books, 
we refer the reader to, say for details of the varia- 
tional method. We first determine the energy functional. 



Ev = {i'trix) I Had I i'tr{x)) 
b kl ^ 



(33) 



4 ■ 4& ■ V2& t^(-l/2,0,6fc|/2a) (34) 



This can be exactly solved for b but leads to rather com- 
plicated terms. As an aside, we also point out that us- 
ing asymptotic forms for U{.) and A'o(.) as z — s- 00, i.e, 
a — > 0, we can recover the position eigenstate obtained 
in Eq. 

From this point, we specialise to our case fci = /s2 = 1. 
Then, for a >> 1, the second term in Eq. (|38|1 can be 
neglected since b increases monotonously with increasing 
a. Now, setting dEy/db = 0, we obtain for a >> 1, 



v) 



1/3 



(39) 



In general, b is dependent on fci and k2 as evident from 
the exact solution to Eq. H38() but we have taken ki — 
^2 = 1. Plugging in this b in Eq. (|32|1 . we obtain the 
adiabatic ground state valid for a >> 1 as. 



(40) 



where fix = Vl + 2ax'^ since we have taken fci = fc2 = 1 
to make the analysis tractable. It is correctly normalised 
but does not give the correct limit for a = 0. It does not 
possess the symmetry a; <-> y. 
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1. 



Position entropy for a >> 1 



The probability density in position representation is 



(41) 



From this, the information entropy with ki — k2 — I can 
be evaluated (see Eq. I|23|) as, 



Sg{a) = l + lnTT-lnVfo- 



I: 



(42) 



where I2 is given by. 




b f 3 b 



(43) 



and erfi(.) is the imaginary error function and 2-^2 (•) is 
the generalised hypergeometric function. In this form, 
the information entropy in Eq. 142(1 is not particularly 
illuminating. For a » 1, the contributions from 2nd 
and 4th term in Eq. (|43|l is negligible. Thus, we obtain 
a simplified form for the spatial entropy, 

4"4^n2fc 
Now, substituting for b from Eq. H38|l . we obtain 



^,(a) = l + ln^-lnV6+^--ln(-) (44) 



Sq{a) = 1 + InTT 



--Ina 



(45) 



Thus, for large a, the ground state entropy falls loga- 
rithmically with a slope 1/3. This is verified numerically 
through exact calculations. In Fig. E^b), we show the po- 
sition entropy plotted as a function of ln(Q!). A linear re- 
gression gives the best fit line as Sq = 2.1368—0.3033 Ina, 
verifying the approximate theoretical slope in Eq. ((45|l . 
The intercept is approximately 1 -)- InTr, the entropy of 
the unperturbed oscillator, but we stress that we can- 
not recover a = result from Eq. H45|) . We also point 
out a systematic difference seen in Fig. 01a,b) between 
the theoretical estimate of entropy (crosses) given by Eq. 
(|45|l and the numerical result (circles). This is due to the 
adiabatic approximation becoming less accurate near the 
origin of the potential and this is showing up in the re- 
sults. 
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FIG. 4: (a) Entropy in momentum representation and in (b) 
spatial representation for the Hamiltonian in Eq. Q with 
fci = ^2 = 1 as a function of a ranging from 30 to 90. Dots 
are the numerical entropies and solid line is the best fit to 
the numerical data. Crosses represent the theoretical result 
based on the adiabatic approach. 



integral in Eq. (|16|l using the ground state in Eq. I|40|) . 
Once again, by transforming variables, we can reduce it 
to a one-dimensional integral, 



51^ 



dx 



1/4 



(46) 



where Q,x is given by Eq. JHJ and b is the variational 
parameter determined in Eq. 1(38(1 . This integral could 
only be performed numerically. Hence, we first determine 
the adiabatic momentum eigenstate by numerically inte- 
grating Eq. 1(46(1 from which the entropies are computed. 
The entropies obtained numerically by the adiabatic ap- 
proach, shown as crosses in Fig^fa), and those obtained 
numerically without any approximation (circles in Fig 
^a)) are in fair agreement with each other for the range 
of a considered in this work. A linear regression gives 
the line of best fit, 



Sp{a) = 2.1754 + 0.3134 Ina 



(47) 



Momentum entropy increases logarithmically with a. We 
notice from our numerical calculations, exact as well as 
the one based on adiabatic approach, that the slope in 
the equation of best fit for the momentum entropy is 
consistently higher than the one for the position entropy. 
This difference accounts for the behaviour of the entropic 
sum presented in the next section. 



C. Entropic sum for large a 



B. Momentum eigenstate and entropy 

We obtain adiabatic ground state in the momentum 
representation for a >> 1 by Fourier transforming the 
position eigenstate in Eq. ((40(1 . We need to perform the 



Finally, we put together the results to look at the 
entropic sum for large a. The Fig [S] shows the en- 
tropic sum for the ground state for a range of parameters 
a — 30, ....90. The best fit line (solid line) is given by, 

S = 4.3249 -I- 0.00708 In a (48) 
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ln(a) 

FIG. 5: (a) Entropy sum for the ground state of the Hamil- 
tonian in Eq. Q with fei = ^2 = 1 as a function of a. Circles 
are numerical values and solid line is the best fit to the nu- 
merical data as given by the Eq. (UHl . 



Firstly, it satisfies the entropic uncertainty inequality in 
Eq. Q). The positive slope in this equation shows that 
the entropic sum increases with the parameter but rather 
slowly. It must also be pointed out that this slope in 
Eq. corresponds closely to the difference in slopes 

corresponding to the best fit equations for position and 
momentum entropy. 

The functional form of Eq. H48(l is quite analogous 
to the entropic sum in Eg. Q conjectured for atoms, 
clusters and nuclei Td, Ti] on the basis of theoretical 
arguments and numerical results. Such logarithmic in- 
crease in entropic sum is noted earlier for ID power-law 
potentials in the semiclassical limit 9]. Thus, this result 
provides an approximate (semi-)thcrcotical approach to 
see the emergence of the functional form in Eq. I|48|) in 
a simple coupled system. 



VI. HIGHLY EXCITED STATES 
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FIG. 6: Spatial (circle) and momentum (square) entropy for 
a highly excited state characterised by approximate quantum 
number pair (110,0) with fci = ^2 = 1 as a function of a. 

then follow the cigcnstate as a function of a. The cal- 
culations are also somewhat cumbersome and hence we 
sample only a small parametric window. 

The results, obtained from exact numerical basis-set 
diagonalisation and then a subsequent entropy calcula- 
tion, shown in the Fig. IHlreveal that the spatial and mo- 
mentum entropies display similar trend as a function of 
a. This is qualitatively different from that of the ground 
state. From the figure, it is also clear that the variation 
with a is not monotonous. The entropic sum evidently 
obeys the inequality in Eq. It is known that for 

Hamiltonian in Eq. © there are certain special peri- 
odic orbits which are responsible for supporting regular 
states in the system j22|. The probability density struc- 
tures are related to the classical phase space structures 
|23| and hence the information entropy is also intimately 
connected to the properties of the local classical phase 
space structures that support such regular states. Thus 
the entropy is modulated by the qualitative nature of 
local classical dynamics as a function of a ^2^. 



In this section, we look at the grey areas that are not 
yet clearly understood. We briefly discuss the EUR for 
the highly excited states. Most eigenstates of chaotic sys- 
tems in the semiclassical limit are irregular or random 
looking states that could be modelled by random matrix 
theory |l9l |. Such random matrix averages represent a 
particular limit at which the system's entropy becomes 
independent of the nonlinearity parameter. Hence we fo- 
cus our attention on the highly excited 'regular' states 
which are characterised by N quanta of excitation in one 
degree of freedom and in the other, often called the 
Born-Oppenheimer type of states 0, Thus they 

could be approximately labelled by the quantum number 
pair {N, 0; a = ao), where >> 1 For these states, 
we will numerically explore the variation of entropies as a 
function of a. In a sense, the entropies of highly excited 
states have already been reported before in a different 
context |23|. Firstly, we choose a particular eigenstate 
characterised by {N,0;ao) , in our case (110,0;0.0), and 



VII. CONCLUSIONS AND DISCUSSIONS 

This work can be broadly divided into two parts. Fo- 
cussing on the ground state of a coupled system, firstly, 
we show, within the framework of adiabatic theory, that 
spatial and momentum eigenstates can be constructed. 
Secondly, using the adiabatic eigenstates we obtain ap- 
proximate results for the entropic uncertainty relations. 
For the ground state, we show that the spatial entropy 
decreases as a function of the coupling parameter a in the 
system while the momentum entropy increases. However, 
the entropic sum increases with the coupling parameter 
a and gets saturated only for alpha = 0. This is reminis- 
cent of the results reported numerically for the entropic 
sums of atoms, molecules and clusters Thus, this 

work provides an approach to obtain analytical results 
for the information entropies and their sums, for coupled 
nonlinear systems as a function of coupling parameter. 



8 



However, for highly excited states, the spatial and mo- 
mentum entropy qualitatively look similar and seem to be 
related to the qualitative nature of local classical dynam- 
ics. In general, the relation between the classical phase 
space structures and the quantum eigenstates is not com- 
pletely understood yet. In this context, one might men- 
tion the semiclassical approaches due to Berry and Bo- 
gomolny 23] based on the Gutzwiller's Trace Formula. 
Thus, we only point out the intricacies involved in the en- 
tropic sums for the highly excited states. In due course, 
it might become possible to relate information entropy 
to classical orbits through these approaches. There are 
already empirical results pointing to this connection |2Clj| . 

The adiabatic theory based approach presented here 
is not without demerits. Strictly speaking, the adiabatic 
effects take over when there is clear separation in the 



time-scales of motion in the two modes that constitute 
the system. Thus, in the case of Hamiltonian in Eq. ((HJ, 
the potential develops channels (see Fig. ^ for large 
values of a and thus facilitates adiabatic effects. The 
adiabatic theory is particularly suitable if the probabil- 
ity density develops structures within these channels, as 
it happens for some of the highly excited states. Since 
the ground state does not enter the channel for large a, 
it is likely to become less accurate for very large values of 
a. Further, this work leaves an interesting question unan- 
swered. What happens in the limit a = oo ? The present 
method of treatment does not seem adequate to answer 
this question. We hope the results here could stimulate 
research on more rigorous approaches to entropic sums 
in a wide variety of coupled systems. 
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